Relation between flux formation and pairing in doped antiferromagnets 
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We demonstrate that patterns formed by the current-current correlation function are landmarks 
which indicate that spin bipolarons form in doped antiferromagnets. Holes which constitute a spin 
bipolaron reside at opposite ends of a line (string) formed by the defects in the antiferromagnetic 
spin background. The string is relatively highly mobile, because the motion of a hole at its end does 
not raise extensively the number of defects, provided that the hole at the other end of the line follows 
along the same track. Appropriate coherent combinations of string states realize some irreducible 
representations of the point group dv Creep of strings favors d- and p-wave states. Some more 
subtle processes decide the symmetry of pairing. The pattern of the current correlation function, 
that defines the structure of fiux, emerges from motion of holes at string ends and coherence factors 
with which string states appear in the wave function of the bound state. Condensation of bipolarons 
and phase coherence between them puts to infinity the correlation length of the current correlation 
function and establishes the flux in the system. 



PACS numbers: 74.25. Jb, 71.10.Fd, 71.27.+E 



I. INTRODUCTION 

Flux phase appeared in the mean-field approaches 
to undoped aad doped antiferromagnets described by 
the t-J modelEJ. The physical meaning of that phase 
and its relation, with pairing was not clear. In some 
announcementscTLl the competition or coexistence be- 
tween superconductivity and the flux phase have been 
discussed in the framework of the mean filed approach 
and the 1/A^ expansion. Some researchers tend to under- 



stand the spin gap phase. Jn doped antiferromagnets as 
the staggered flux phasmla. Other believe that pairing 
grows out of the 7r-flux phaseO or that pew:^ are formed 
by holes circulating in opposite directionst2l. In this pa- 
per we will present a different point of view that the pat- 
tern of currents and pairing originate in formation of spin 
bipolarons of a particular form. We will give some math- 
ematical arguments, based on understanding of pairing in 
weakly doped antiferromagnets, iiLsupport of that idea. 

Recently, Ivanov, Lee and WenE£l reported staggered- 
vorticity correlations of the current-current correlation 
function (jijjki) (CCCF) in the d-wave variational func- 
tions for a weakly doped antiferromagnet, where jij de- 
notes the current flowing on the bond {i,j). Similar pat- 
tern of the CCCF was found for a bound state of tatp 
holes in an exact diagonalization performed by Leungllil. 
Since the mechanism of hole binding in weakly doped 
antiferromagnets is relatively well understoodtj, it is a 
tempting task to resolve whether there exist a deep rela- 
tion between binding and the formation of flux patterns. 

Fast motion, with the rate ~ i, of a hole created at 
some site in the Neel state—forms defects in the origi- 
nal arrangement of spinaljo. These defects, that are 
spins which have been turned upside down and may also 
be called magnons, lie on the track of the hole. Strings 



pinned to a central site form a potential well for a hole 
which conflnes its motion in the vicinity of this site. 
Much slower processes related to the inversion of anti- 
parallel spins at nearest neighbor sites which occur at 
the rate ~ J shorten strings and give rise to coherent 
propagation of the hole. Since the separation of energy 
scales for hopping of a hole and annihilation of defects 
(magnons) is pronounced it is plausible to introduce the 
notion of a spin polaron. We may define in this context 
a spin polaron as a solution of the, Schrodinger equation 
for a particle in the potential welt3. Slower processes are 
neglected at this stage of analysis. "Orbital" states of po- 
larons created at all possible sites exhaust, in principle, 
the relevant portion of the Hilbert space. In practice, the 
calculation may be confined to some low excited states 
or even to the polaron groundstate. By calculating the 
matrix elements of the full Hamiltonian in the polaron 
basis we derive an effective Hamiltonian expressed in the 
polaron language. Some processes, like oscillations of the 
hole in the vicinity of a polaron center, are already incor- 
porated into the eigenenergy of polaronic states. The rest 
either renormalizes the eigenenergies or gives rise to off- 
diagonal matrix elements in the polaronic Hamiltonian. 
The latter circumstance occurs in the case of shortening 
of strings by an appropriate term in the original Hamil- 
tonian. The approach outlined above has some obvious 
limitations. It is applicable provided that the correlation 
length of background antiferromagnetic spin correlations 
is larger than the polaron radius. 

We extend now the concept of spin-^larons to the 
description of interaction between holcsll3. Since we are 
going to apply this construction to the calculation of the 
current correlation function, we will provide the reader 
with more details. We define a spin bipolaron j^I'^^ j)) 
as a combination of states which may be reached by in- 



dependent hopping of holes created at a pair of nearest 
neighbor sites (i,j}, 






(1) 



Vi parameterizes the geometry of a path along which the 
hole has been moving and IPijVj) is a state which has 
been created in this way. Due to Fermi statistics, order 
in which holes have been originally created in {"^(^ijj) is 
relevant. We adopt a convention, according to which a 
hole is created first at the site i and i belongs to the 
even sublattice. At this stage of considerations we pro- 
hibit, by definition, each hole to follow along the trace 
left by the accompanying hole. By means of that some- 
what artificial restriction we achieve that spin bipolarons 
are localized and we may proceed as in the case of a sin- 
gle hole. We make a further approximation and neglect 
all path details including the possibility of path crossing. 
It is self-evident that the coefficients a-p^^-p. for a bipo- 
laron in the groundstate will depend only on the lengths 
/i, ly of paths Vi and Vj, a-Pi,Vj = a/^,iy. The Schrodinger 
equation that describes a hole pair in a potential well and 
defines the spin bipolaron, may be written as 



t [ctf^-i.u + {z - l)ta^+i,^ -I- a^,i._i + (z - l)to^,^ 



+1J 



+J I 4 + /i + i/- -(5^,^ 1 a^^^ == £;2a^,i., (2) 

where a^^i, — for ^ <Q oi v < and z=4. The form of 
this equation is easy to understand. The first expression 
in it is related to the fact that each path may be reached 
from one shorter path and (z-1) different longer paths by 
a hop of the hole. The second term counts the number of 
pairs of nearest neighbor sites, that are not occupied by 
anti-parallel spins. Every such a "broken bond" raises 
the energy by J/2 in comparison to the energy of the 
Neel state. 

Formation of bipolarons, and, in particular, strings 
that connect holes, is an effective way of lowering energy. 
A compromise between two opposite tendencies to min- 
imize the kinetic energy of holes and to reduce the dis- 
turbance of the antiferromagnetic background is reached 
by means of that process. By shrinking at one end and 
expanding at the opposite-end, a string may move, while 
keeping a moderate lengthliH. The application of the spin- 
polaron scenario to binding of holes in doped antiferro- 
magnets gave rise to better than qualitatpe agreement 
with results of numerical diagonalizationsllS. In partic- 
ular the hierarchy of states that realize irreducible rep- 
resentations for all momenta and symmetries allowed Jay 
the geometry of the 4x4 cluster have been reproducedllj. 
It turns out that creep of strings favors d-wave and p- 
wave states. More subtle processes decide the symmetry 
of the bound state in favor of the d-wave. Some modifica- 
tions of the t-J model, like asymmetry of the magaatic in- 
teraction may lead to stability of the p-wave stateEj. The 
weight of bipolaronic states in the bound state amounts 



to about 80% in tiie case of the d-wave and 50% in the 
case of the p-wavell3. Thus it is apparent that in a qual- 
itative description of current correlation function in the 
presence of pairing we may neglect monopolaronic states. 
S-, d- and p-wave symmetries piay be realized as co- 
herent sums of bipolaronic statesES: ^u j\ ^(i,j)\^ (i,j))i 
where Suj\ = 1 if (j,j) is horizontal and Suj\ = —1 
provided that {i,i) is vertical for the d-wave, while for 
the p-wave Si^i^j) vanishes if {i,j) is vertical, St^ij) = 1 
provided that j is on the right side of i and Si^i^j) = — 1 
if on the left. In order to make the story short we have 
skipped many details. The full construction of bipolarons 
that conform with irreducible representationSpOf the point 
group C4„ may be found in an earlier articleta. 



II. FLUX PATTERN OF THE D-WAVE PAIRED 
STATE 

As a first example of a contribution to the CCCF we 
consider a process shown in Fig. |l|(a). The middle row 
and the lowest row depict string states of the type \Vi, Vj) 
which are coupled by action of a product of current op- 
erators on outer bonds. Circles represents positions of 
bipolarons, or in other words sites at which holes have 
been created, while thin arrows symbolize paths. In the 
case of the state in the second row, a single defect in 
the Neel structure occupies the second site from the left, 
while in the case of the final state in the third row a 
single defect (magnon) occupies the third site from the 
left. Current operators which are given by the formula 
jij oc i{c\^Cja — Cju'^irr) shift holcs in the middle row. 
Holes are moved in the same direction, which is repre- 
sented by the arrows in the uppermost row. 
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FIG. 1. Some processes which contribute to the CCCF 
function when holes are shifted by a product of current oper- 
ators in the same direction on outer bonds. 

If we adopt a convention that the direction of current 
on both outer bonds is from left to right the contribution 
to the CCCF is t^a^ i- The factor i"^ which stems from 
product of two current operators has been absorbed by 
a change of sign related to the fact that the hole cre- 
ated at the even site in the middle row landed at a site 
which is occupied in the lowest row by a hole created at 
the odd site. Fig. Wa) is an example of a contribution 
to a the CCCF for two bonds which may be connected 
by a line consisting of a single bond. If we incorporate 
the outer bonds to that line the length of it amounts to 
three. Since we expect that the absolute value of am- 
plitudes Qf^^jy decreases when the total length of two in- 
volved paths ^ + V increases, only the shortest possible 



string states will decide the sign of the CCCF for two 
chosen bonds. Figures || (b), (c), (d) depict the rest of 
processes which involve shortest paths for holes which are 
moved from left to right at outer bonds of a line which 
length is 3. The contributions from these processes to 
the CCCF are — i^a§j, —t^a^ i, and t^a^ i respectively. 
In the final formula these values are multiplied by factors 
Siij\ with which bipolaronic states at points represented 
by open and solid circles appear in the coherent sum that 
defines the bound state. Analogous processes contribute 
to the CCCF if holes are shifted from right to left on 
both outer bonds. 

Fig. depicts contributions to the CCCF for holes 
which move inward on outer bonds of a line of length 
3. 
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FIG. 2. Some processes which contribute to the CCCF 
when holes are moved by a product of current operators in 
opposite directions on outer bonds. 

In the original string state, in the second row of Fig. g 
(a), to which the product of current operators is applied, 
holes occupy outer sites, while two magnons reside on 
a pair of sites in the middle. These magnons (defects) 
disappear when holes are shifted to the central position 
in the lowest row. Since holes that have been created in 
even and odd sublattices, exchange positions, the Fermi 
statistics again contributes the factor "—1" but holes are 
moved in opposite directions, and the contribution to the 
CCCF becomes — i^ao.oao,2- Amplitudes of processes in 
Figures g (b) and (c) are t^ao,oQ!i,i and —i^ao, 0^0,2 re- 
spectively. The directions of hole motion may be inverted 
if the product of current operators is applied to the string 
state in the lowest row which gives rise to next three con- 
tributions to the CCCF on outer bonds. We can use the 
insights gained from the analysis of a line which consists 
of 3 links in total, to obtain a formula for the contribution 
from bipolaronic states which lie on a line of arbitrary 
length Z, to the CCCF on outer bonds of that line, 

2 Y^ X^ ( _^\m+n+l 



\~^) Jm^nt Oim-l,l-m-lOin-2,l- 



m— 1 n—2 

+ 2 2J /X~'^)™'^"'SmSnt'^am-l,l-7nan-2,l-n-l, (3) 
m — 1 n—2 

where Sm is a coherence factor with which the bipolaron 
at m-th bond in the line appears in a coherent sum that 
represents a bound state. An additional lesson of the 
above analysis is that the shortest string states which 
contribute with the highest weight to the CCCF for a 
pair of bonds, lie on shortest paths that lead through 



lattice points and incorporate these bonds. If the length 
of such paths amounts to I, the length of involved strings 
is / — 1, when holes are moved by a product of current 
operators in the same direction or I and 1 — 2, when the 
holes are shifted in opposite directions. As we have al- 
ready mentioned, the leading contribution to the CCCF 
for a pair of chosen bonds is the sum of the expression 
(^ over all shortest paths that incorporate those bonds. 
Fig. represents in units 10~^ the CCCF divided by the 
hole concentration, that has been calculated in this way 
for the d-wave bound state of two holes. Arrows rep- 
resent the direction of current on respective bonds. We 
notice that despite the crudeness of the applied approx- 
imation, agreement sath results of exact diagonalization 
performed by LeungEll is much better than qualitative. It 
is an intriguing question, if that coincidence is a manifes- 
tation of a deeper relation between pairing and patterns 
of the CCCF. In order to see that this is the case we 
assume that the logarithm of the weight of a string a 
is proportional to its length I. That assumption is very 
natural because by expanding the string by one lattice 
spacing we raise the number of broken bonds by 2 and 
reach 6 longer string states, independently of the value 
of I. In our description of spin polarons parameters cu^^^, 
where /i + i^ = / — 1, are related with strings of length I, 
but in general a^,i/ ^ a.f_i' y for ^, + v = fi' + v', which 
is the prize which we pay for pinning strings and making 
bipolarons localized. 
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FIG. 3. {jkijm.n)/x in units 10""", where x is the hole 
concentration, at J = 0.3t. The reference bond has been 
marked by a circle. 

We notice that the total length of strings contributing 
to elements of sums in (H) is 21 — 2. Under the previ- 
ously made assumption, that q;^,,^ oc exp(C(/i + i^)), all 
products of two weights a are the same in both sums. If 
we neglect them and evaluate the sum of expressions for 
all shortest lines that link bonds we again get the stag- 
gered flux pattern presented in Fig. H. We notice that 
the CCCF defined in this way has the expected structure 
and apart from bonds attached to the reference bond the 
CCCF is given by twice the number of shortest paths 



that connect bonds with the reference bond. This is not 
by coincidence, because the sum (|3|) with factors a ne- 
glected, may be evaluated explicitly, 



i-\ I 



2 Y. ^(-l)™+"+i5™5„ 

771—1 n— 2 
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This formula implies that the pattern of the CCCF is de- 
termined in a simple way by the symmetry of the bound 
state or to be more specific by the coherence factors with 
which bipolarons on involved bonds appear in the wave 
function. 
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FIG. 4. The sum of expressions (|3|) for all shortest paths 
that connect bonds with a bond marked by an open arrow for 
the d-wave symmetry. The weights a have been neglected in 
the evaluation. 



DDCF obtained for the d-wave state by adding contri- 
butions (g) from all shortest lines that connect relevant 
sites. 
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FIG. 5. Contributions to the DDCF, Chh{r) at J = 0.3t 
from string states that he on shortest paths which connect 
relevant sites. 

Notwithstanding the simplicity of our approach the 
similarity with results of numerical diagonalization by 
Leungllll is evident. Oscillations of Chhir) with distance 
and minima at points which lie on an axis are visible in 
both calculations. Thus, we draw a conclusion that the 
pattern of the CCCF, the structure of the DDCF and the 
relation between rates of their decay with distance in the 
rf-wave state is related with formation of bipolarons and 
strings. 



III. DENSITY CORRELATION FUNCTION 

There exists another quantity which behavior is de- 
termined by existence of spin bipolarons. Density 
matrix renormaiiaation group calculationsE3 and exact 
diagonalizationsllil point out that the density-density cor- 
relation function (DDCF) Chhir) for holes as a function 
of distance r between them reveals a characteristic struc- 
ture in the d-wave bound state. The most pronounced 
feature of it is that the DDCF decays much faster along 
lines in the x (y) direction, especially if the distance r 
is an even number. We assume again that the highest 
contributions to the DDCF stem from shortest possible 
string states. These strings should stretch between rele- 
vant sites in order to contribute to the DDCF. Thus, the 
decay of the DDCF with distance is related to the decay 
of CCCF by the dependence of weights parameters a on 
string length. Gray disks in Fig. |j point out to the rela- 
tion between the DDCF at a pair of sites and CCCF on 
bonds attached to them. After a short evaluation, the 
contribution to the DDCF from strings that lie on a line 
takes the form. 
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IV. FLUX PATTERN OF THE P-WAVE PAIRED 
STATE 

Motivated by this observation we may try to guess the 
pattern of CCCF in the p-wave bound state. A state 
with that symmetry may became a groundstatc if the 
exchange interaction between spins in the t-J model be- 
comes anisotropic. By applying the formula (0) we get 
the structure presented in Fig. || for the case of the sym- 
metry p^. 



For the d-wave state that sum exactly vanishes for points 
at one of the axis (x or y) if / is even. Fig. ^ depicts the 



FIG. 6. The pattern of currents in the CCCF for the 
Pi^-wave symmetry. 

No current flows in the vertical direction because there 
are no bipolarons on vertical bonds. The cancellation 
which lead to the formula (0) is not exact in the case 
of the expression (^, and in addition some higher order 
contribution to the CCCF from spin bipolarons will ap- 
pear for longer paths. Evaluation of these contributions 



to the CCCF is beyond the scope of this paper. Monopo- 
larons may also obscure the picture presented in Fig. g 
because the weight of bipolarons in the p-wave state is 
smaller than for the d-wave state. 



V. CONCLUSIONS 

In summary, we have demonstrated that the vortic- 
ity correlations in weakly doped antiferromagnets have 
their origin in formation of spin bipolarons and strings. 
The flux pattern is determined by the symmetry of the 
bound state. The rate of decay with distance of differ- 
ent correlation functions related to current or density is 
determined by the rate with which the weight of strings 
diminishes with their length. That observation indicates 
that pairing and vorticity correlations in doped antiferro- 
magnets have common origin in formation of bipolarons 
or in other words of strings built by magnons that con- 
nect holes at their endpoints. 

The proximity of both states war already observed in 
mean field calculations. Conclusions drawn from that 
observation that pairing grows out of flux or circulating 
currents are questionable in the light casted by calcula- 
tions presented in previous sections. Also the statement 
that pairs are formed by holes that circulate in opposite 
directions does not agree with our scenario. It seems that 
agreement between analytical calculations and exact di- 
agonalizations suggest that holes that form a bound state 
oscillate in all possible directions at opposite endpoints 
of a string. 

Depletion of the single particle spectral weight around 
the Fermi levelEj known as pseudogap and exceptional 
non-FeLmi^liquid (NFL)properties of copper-oxide based 
metalst3'E2l are characteristic features of those systems 
that decide about the complexity of their physics. Some 
suggestions have been expressed that phase fiuctuations 
are responsible for the pseudogap phenomenonnJ. These 
ideas are tightly rekd:ed with the Bose-Einstein conden- 
sation of local pair£3. The microscopic derivation in the 
framework of a two baiid model of the marginal Fermi 
liquid (MFL) conceptE3 which seems to describe quite 
well the NFL behavior observed in copper oxides is based 
on existence of circulating currents. The calculation pre- 
sented in this paper suggests that both local pairs respon- 
sible for pseudogap phenomena and circulating currents 
which give rise to the MFL behavior may be attributed to 
formation of spin bipolarons. Finally, it is important to 
emphasize that both phenomena are tightly related and 
can not be discussed separately in the case of the copper 
oxides, which is also in agreement with the suggestion 
their relation with spin bipolarons and string states. 
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